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I .  INTRODUCTION 


The  simplest  theory  for  nonlinear  wave  propagation  in  a  straight, 
uniform  rod  may  be  written  as  follows. 
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(1) 


The  axial  stress,  S,  depends  only  on  the  axial  strain,  e,  the  axial 
particle  velocity  is  v,  and  the  reference  density,  p,  is  a  constant.  The 
axial  coordinate  along  the  rod  is  Z,  t  is  time,  and  subscripts  denote 
partial  differentiation.  In  this  version  of  rod  theory,  which  is 
formally  the  same  as  one  dimensional,  inviscid  gas  dynamics,  the  char¬ 
acteristic  speed  of  propagation,  c,  is  given  by  pc2  «  dS/de,  and  dis¬ 
continuities  in  stress,  strain,  or  velocity  all  travel  at  the  shock 
speed  given  by  pc2  =  [S]/[e],  where  [.]  denotes  the  jump  in  a  quantity 
across  the  shock. 


Equation  (1)  has  been  used  to  interpret  wave  propagation  experiments 
in  which  the  strain-time  profile  is  measured  at  several  stations  in  the 
r.od.1*^  jf  each  level  of  strain  propagates  at  a  constant  velocity, 
c  (e) ,  then  both  stress  and  particle  velocity  may  be  computed  as  a  func¬ 
tion  of  strain  as  follows. 


fe(Z,t) 

v  =  v  -  I  c  de 

0  Je  P 

o 

(2) 

«e(Z,t)  2 

S  =  S  +  1  pc  de 

0  Je  P 


The  initial  values  ahead  of  the  wave  are  vQ,  SQ  and  e^.  Equation  (2)  is 
a  simple  wave  solution  to  (1) ,  exhibiting  amplitude  dispersion,  but  not 
geometric  dispersion. 


\j.  F.  Bell,  The  Experimental  Foundations  of  Solid  Mechanics,  Handbuch 
der  Physik,  Vol.  VIa/1,  Springer- Ver lag.  New  York,  1973. 

2 

G.  £.  Hauver,  "Penetration  with  Instrumented  Rods,"  Int.  J.  Eng.  Sci, 
16,  1978,  pp.  871-877. 
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A  rod,  no  matter  how  slender,  is  actually  a  three-dimensional 
object,  of  course.  The  dispersive  influence  of  a  finite  diameter  is  well 
known  in  linear  theories,  eg.3*4  It  is  less  well  known  in  nonlinear 
theories5'6  where  it  may  be  expected  to  play  a  role  as  well,  whenever 
the  length  scale  in  a  wave  pulse  is  comparable  in  magnitude  to  the  rod 
diameter.  In  this  paper  the  effects  of  finite  lateral  dimensions  will 
be  considered  by  modeling  a  rod  as  a  one-dimensional  elastic  structure 
with  one  internal  variable  used  to  represent  the  transverse,  axisymmetric 
motion.  Such  a  structure  is  a  special  case  of  an  intrinsic  rod  theory 
as  described  by  Antman7. 


II.  ONE-DIMENSIONAL  EQUATIONS  AND  RELATIONSHIP  TO 
THREE-DIMENSIONAL  ELASTICITY 

A  straight  cylindrical  rod  of  radius  a  is  assumed  to  have  an  elastic 
stored  energy  density  per  unit  length 

W  -  ir  a2  W(w*  ,u,u') .  (3) 

The  axial  displacement  is  w,  u  is  a  measure  of  radial  strain,  and  the 
dash  represents  differentiation  with  respect  to  the  axial  coordinate,  Z. 
The  kinetic  energy  density  per  unit  length  is  given  by 

K  -  *  a2^*2  +  *jp2u2),  (4) 

where  and  p2  are  the  appropriate  mass  densities  (made  definite  by 
equation  (13))  for  the  axial  and  radial  motions  respectively,  and  the  dot 
denotes  differentiation  with  respect  to  time.  The  Euler-LaGrange  equa¬ 
tions  corresponding  to  (3)  and  (4)  are  given  by  (5) . 


3R.  Skalak,  "Longitudinal  Impact  of  a  Semi-Infinite  Circular  Elastic  Bar," 
J.  Appl.  Mech.,  24,  1957,  pp.  59-64. 

4W.  A.  Green,  "Dispersion  Relations  for  Elastic  Waves  in  Bars,"  in 
Progress  in  Solid  Mechanics,  Vol.  I,  ed.  I.  N.  Sneddon  and  R.  Hill, 
North-Holland,  Amsterdam,  1960. 

5J.  H.  Shea,  "Propagation  of  Plastic  Strain  Pulses  in  Cylindrical  Lead 
Bars,"  J.  Appl,  Phys.,  39,  1968,  pp.  4004-4011. 

6G.  P.  DeVault,  "The  Effect  of  Lateral  Inertia  on  the  Propagation  of 
Plastic  Strain  in  a  Cylindrical  Rod,"  J.  Mech.  Phys.  Sol.,  13,  1965, 
pp.  55-68. 

7S.  S.  Aritman,  The  Theory  of  Rods,  Hanbuch  der  Physik,  Vol.  VIa/2, 
Springer-Verlag,  New  York,  1972. 
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s'  -  P]LV» 

(5 

Q'  -  P  -  p2u 

The  forces  S,  P,  and  Q  are  obtained  from  the  stored  energy  as  follows. 
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aw  _  aw 

aw*"  *  p  “  3u 


(6) 


Equations  (5)  and  (6)  may  be  interpreted  in  terms  of  the  three- 
dimensional  theory  of  nonlinear  elasticity.  For  the  axisymmetric  motions 
considered  here,  the  deformation  from  an  unstressed  reference  configura¬ 
tion  (R,0, Z)  to  the  present  configuration  (r,9,z)  may  be  written 


r  =  r(R, Z,t) ,  0  =  0,  z  «  z(R,Z,t)  ,  (7) 

and  the  corresponding  equations  of  motion  are  given  by 
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(8) 


The  stress  components  used  here  are  obtained  from  the  unsymmetric  Piola- 
Kirchhoff  tensor  referred  to  unit  basis  vectors  in  the  given  cylindrical 
coordinates  (r,0,z)  and  (R,0,Z). 

T  =  T1Q  gj  S  Ea  (9) 

If  (8)  is  multiplied  by  R,  and  if  the  resulting  equation  and  (8)2  are 
both  averaged  over  the  cross-section,  equations  (10)  are  obtained  for  the 
case  of  stress-free  lateral  surfaces. 


-4-  si  pzl « 

ira 

w  dA-^f ^ 


+  T00)dA 


(10) 


Now  if  the  following  identifications  are  made 
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ira 

f RTrZ  dA  , 

and  if  a  first  approximation  for  the  position  functions  r  and  z  is 
assumed  to  be 


r  =  R[1  +  u(Z,t)]  , 
z  =  Z  +  w(Z,t)  , 


(12) 


then  equations  (5)  are  obtained  with 

2 

P1  =  p  '  p2  =  iip& 


(13) 


In  all  of  the  preceding  discussion  it  has  been  tacitly  assumed  that 
the  Z-axis  is  an  axis  of  material  symmetry  so  that  no  angular  motion  will 
occur  anywhere  in  the  cross-section.  Furthermore,  since  the  strain 
energy  must  be  invariant  under  reversal  of  the  z-axis,  and  if  it  is 
assumed  to  be  invariant  under  reversal  of  the  Z-axis  as  well,  W  must  have 
the  property 


W(w',u,u')  =  W(w',u,-u') 


(14) 


It  follows  that  Q  is  odd  in  u'  and  S  and  P  are  even  in  u' . 

It  should  be  remarked  that  equations  (5)  and  (6)  have  been  used  to 
describe  one-dimensional  waves  in  porous  materials®*9* 10,  and  it  has 


fl 

J.  W.  Nunziato  and  E.  K.  Walsh,  "On  the  Influence  of  Void  Compaction  and 
Material  Non-uniformity  on  the  Propagation  of  One-Dimensional  Accelera¬ 
tion  Waves  in  Granular  Materials,"  Arch.  Rat.  Mech.  Anal.,  64,  1977, 
pp.  299-316  and  Adendum,  Arch.  Rat.  Mech.  Anal.,  67,  1977,  pp.  395-397. 

q 

J.  W.  Nunziato  and  E.  K.  Walsh,  "One-Dimensional  Shock  Waves  in  Uni¬ 
formly  Distributed  Granular  Materials,"  Int.  J.  Solids  and  Structures, 
14,  1978,  pp.  681-689. 

10D.  F.  Parker  and  B.  R.  Seymour,  "Finite  Amplitude  One-Dimensional 
Pulses  in  an  Inhomogeneous  Granular  Material,"  Arch.  Rat.  Mech.  Anal., 
72,  1980,  pp.  265-284. 
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also  been  suggested  that  perhaps  they  could  be  used  to  describe  waves  in 
a  layered  composite  .  There  are  undoubtedly  other  applications  as  well. 


III.  STATIC  SOLUTIONS 

It  is  convenient  to  distinguish  two  types  of  static  solutions.  In 
the  first  type  u'  is  identically  zero,  so  equations  (5)  reduce  to 


S (w ' ,u, 0)  ■  SQ  , 
P(w',u,0)  =  0  . 


(IS) 


Since  the  reference  configuration  has  been  assumed  to  be  unstressed  ,  it 
is  required  that  when  w'  **  u  =  0,  (15)  is  satisfied  with  S  =  0.  It  will 
be  further  assumed  that  the  slope  along  the  curve  defined  by  (15) 2  is 
always  negative  and  bounded. 


du  PW 


v(w' ,u,0)>  0 


(16) 


where  the  subscripts  u  and  w'  denote  partial  differentiation  with  respect 
to  the  arguments  of  P.  This  corresponds  to  the  physically  reasonable 
assumption  that  the  rod  contracts  (expands)  laterally  as  it  is  stretched 
(compressed).  Finally,  it  will  be  assumed  that  S  increases  monotonically 
as  w'  increases  along  the  curve  P(w',u,0)  =  0.  That  is 


Sw'  -  vSu  =  E(w’’u’°)  >  °- 


(17) 


The  bar  modulus,  E,  and  Poisson's  ratio,  v,  will  be  considered  as  defined 

by  (17)  and  (16)  with  the  inequalities  holding  as  well  for  all  values  of 

their  arguments,  not  just  along  the  curve  of  (15)2-  Since  W  has  a  minimum 

at  the  unstressed  reference  configuration  (w'  =  u  =  u'  =  0) ,  Pu  >  0 

there,  so  by  (16)  P  ,  >  0  and  P  >  0  everywhere.  This  also  implies  from 
w  u 

(17)  that  W  ,  ,W  -  W  .  >  0  so  that  necking  instabilities  have  been 
w'w'uu  uw'  6 

ruled  out  (see  Antman  [12],  pp.  97). 

In  the  second  type  of  static  solution  u'  is  not  identically  zero. 
There  are  two  integrals  of  (5) . 


^M.  F.  McCarthy,  private  communication. 

12 

S.  S.  Antman,  "Qualitative  Theory  of  the  Ordinary  Differential  Equa¬ 
tions  of  Nonlinear  Elasticity,"  in  Mechanics  Today,  V.  1,  ed.  S.  Nemat- 
Nasser,  Pergamon,  New  York,  1972. 
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(18) 


S(w',u,u')  -  SQ 

w'QCw'.u.u’)  +  w'So-W(w' ,u,u')  »  B 

If  these  may  be  solved  for  u*  and  w'  as  functions  of  u  and  the  constants 
S9  and  B,  then  u  may  be  found  by  quadrature,  and  the  axial  strain  w'  is 
given  parametrically  through  u.  The  study  of  equations  (18)  will  not  be 
pursued  further  here  except  to  remark  that  the  version  arising  from 
linear  elasticity  gives  boundary  layer  solutions  with  Q  decaying  exponen¬ 
tially  from  the  ends.  Phase  plane  analysis  indicates  that  the  behavior 
of  the  fully  nonlinear  equations  is  similar. 


IV.  SOLUTIONS  OF  THE  LINEARIZED  EQUATIONS 

Before  considering  the  dynamic  nonlinear  equations,  it  will  be  use¬ 
ful  to  examine  some  of  the  properties  of  solutions  of  the  linearized 
equations.  The  three-dimensional  strain  energy  for  an  isotropic  mater¬ 
ial  with  displacement  field  given  by  (12),  when  averaged  over  the  cross- 
section,  leads  to  (19)  for  the  forces. 


S  =  (A+2y)w'  +  2Xu  , 

P  =  4(A+y)u  +  2Aw'  ,  (19) 

r  i-2  1 

Q  *  Ha  yu  , 


where  A  and  y  are  the  usual  Lame  constants.  With  nondimensional  variables 
5  =  Z/a  and  t  =  ct/a,  where  c  is  a  speed  that  is  unspecified  for  now,  and 
with  w  *  w/a,  equations  (5)  become 


W«  +  2  A+2  u; 


u  -  [*8  u  +  4  —  w  1 
«  l  V  H  CJ 


— y  w  , 

.  2  TT 
'1 

2 

V  u 

,2  TT 


(20) 


where  subscripts  again  denote  partial  differentiation.  Subscripts 
indicate  partial  differentiation,  and  the  overbar  has  been  dropped  from 
9.  The  longitudinal  and  shear  wave  speeds  are  c^  ■  / (A+2y) /p  and 

c,  •  /y/p  .  These  equations  in  dimensional  form  were  first  given  by 
Mind 1 in  and  Herrmann13. 


13R.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compres- 
sional  Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S. 
National  Congress  of  Applied  Mechanics,  1950,  pp.  187-191. 
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Since  all  the  coefficients  are  constants,  either  w  or  u  may  be 
eliminated  from  (20)  to  obtain  (21) . 


(21) 


(w,u) 


0 


The  one-dimensional  bar  speed  is  c,  =  /E/p  where  E  is  Young's  modulus, 

E  =  y(3A+  2y)/(X+y).  Equation  (21j  exhibits  what  Whitham  has  called  a 
hierarchy  of  wave  speeds^*  Whereas  discontinuities  can  only  propagate 
with  the  speeds  c^/c  or  C2/C,  it  is  well  known  that  in  a  thin  bar  the 
main  disturbance  in  a  pulse  travels  with  speed  ch/c.  Whitham  has 

that 

the  case  where  the  orders  differ  by  two,  as  here,  also  exhibits  a  stable 
hierarchical  wave  structure,  so  long  as  c2<c^)<c1,  which  is  always  the  case 
for  linear  elastic  rods.  “  '  J 


discussed  extensively  the  hierarchical  case  where  the  orders  of  the 
highest  and  lowest  derivatives  differ  by  one.  Wu*°  has  pointed  out 


In  Whitham's  case  the  higher  order  waves  decay  exponentially,  and 
the  lower  order  wave,  which  carries  the  main  pulse,  has  a  diffusing  front. 
In  the  present  case,  discontinuities  in  the  higher  order  waves  do  not 
decay  with  time,  but  the  pulse  decays  rapidly  behind  the  front.  The  main 
pulse  travels  with  the  lower  order  speed,  and  it  too  has  a  diffusing 
front . 


To  give  substance  to  the  preceding  remarks,  consider  a  boundary  - 
initial  value  problem  for  (20)  . 


B.V.:  u^(o,t)  =0,  w  (0  ,t)  =  f (r) 
I.V. :  u(C ,0)  =  w(C ,0)  =  0 


(22) 


A  straightforward  application  of  the  LaPlace  transform  on  t  leads  to  the 
following  representations  for  u  and  w. 


14 

G.  B.  Whitham,  "Some  Comments  on  Wave  Propagation  and  Shock  Wave  Struc¬ 
ture  with  Applications  to  Magnetohydrodynamics,"  Comm,  on  Pure  and 
Appl .  Math,  12,  1959,  pp.  113-158. 

^G.  B.  Whitham,  Linear  and  Nonlinear  Waves,  John  Wiley,  New  York,  1974. 

16T.  T.  Wu,  "A  Note  on  the  Stability  Condition  for  Certain  Wave  Propaga¬ 
tion  Problems,"  Comm,  on  Pure  and  Appl.  Math.,  14,  1961,  pp.  745-747. 
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1 


?(p)  is  the  transform  of  f(x),  and  m^(p)  andm2(p)  are  the  roots  with 
negative  real  parts  (for  outgoing  waves)  of 


2  2  c 

m-p  -2 

V  c2  / 


0 


(25) 


In  each  equation  of  (23)  the  first  integral  represents  a  wave  whose  front 
travels  at  speed  c1  and  the  second  integral  represents  a  wave  whose  front 
travels  at  speed  c^. 

•  To  examine  the  behavior  of  the  solution  near  the  wavefronts,  let  p 
be  large.  Expansions  for  m^  and  m^  are 


12 


The  leading  terms  in  the  integrals  for  w  may  be  written  asymptotically 
as 


c 

The  first  integral  in  (27)  vanishes  for  x  <  —  c  and  the  second  vanishes 

C1 

for  x  <  —  5.  Thus,  near  the  leading  wavefront  only  the  first  integral 
c2 

is  required,  which  may  be  written,  after  a  little  manipulation  and  use 
of  the  convolution  theorem,  as  follows. 
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JQ  and  are  Bessel  functions.  Suppose  that  f(x)  =  w  0  h(x)  where  h(x) 
is  the  Heavisde  step  function  and  wt°  is  a  constant.  The  change  of 


variable  y  =  r 
ence  17,  pg.  64) 


4 


1  2  c 

T'Z  -  ■£_ 


y  V"  and  use  of  the  identity  (taken  from  Refer- 


Jx(y)  dy 

ST— 2 

Vy  +  n 


now  reduces  (29)  to 


(30) 


Equation  (31)  shows  that  at  points  removed  from  the  boundary  the  velocity 
becomes  exponentially  small  shortly  after  passage  of  the  wave  front. 
Equation  (29)  is  plotted  in  Fig.  1  with  f(x)  =  h(x). 

Similar  analysis  may  be  applied  to  the  other  integrals  in  (23). 

The  second  integral  for  w^  represents  a  wave  whose  front  moves  with 
speed  c2,  and  is  similar  in  form  to  (29),  but  smoothed  by  a  double 
integration  in  x.  u  has  a  contribution  at  each  wavefront  that  is 
similar  to  (29),  but  Smoothed  by  a  single  integration  in  x. 


17F.  Bowman,  "Introduction  to  Bessel  Functions,"  Dover,  New  York,  1958. 
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To  examine  the  main  pulse  in  (23)  it  is  convenient  to  follow 
Whitham's  approach  to  the  problem  (Reference  15,  pg.  345).  First,  the 

exponent  in  one  of  the  integrals  is  written  as  T^p  +  m^p) 

the  saddle  point  method  is  used  to  find  an  asymptotic  expansion  valid 

at  large  t.  For  a  fixed  value  of  C/t  saddles  occur  at 

% 

3p  (p  +  ^CP)  f)  ■  1  +  VCp)  f  =  0.  (32) 

Solution  of  (32)  gives  p  =  p  >  and  the  exponent  is  maximized  with 
respect  to  C/t  when  m^[p(t/T)]  =  0  because  of  (32)  and  (33). 


jjW  [f  *  "l®  I]  ■  -1®  *  [*  *  V®  ?]  I$hr m  0  t33) 

2 

From  (25)  if  m^  =  0,  either  p  =  0  or  p^  *  -8  -j-  •  By  differentiating 

c 

(25)  once  with  respect  to  p,  it  is  clear  that  the  correct  choice  is  p  =  0, 
and  by  differentiating  (25)  a  second  time  with  respect  to  p,  it  is  found 
that 

ml ' (°)  -  ±  §•  »  (34) 


and  so  from  (32) 


(35) 


Equation  (35)  shows  that  the  main  disturbance  travels  with  the  bar 
speed,  c^.  The  correct  asymptotic  expansions  will  be  obtained  by  expand 
ing  the  integrands  in  (23)  about  p  =  0.  The  following  discussion  will 
be  simplified  by  taking  c  =  c^. 

It  has  been  shown  that  m1(0)  ■  0,  but  it  does  not  follow  that 

21  2  4 

n»2(0)  =  0.  Thereforg,  let  m  =  a  +  bp  +  dp  +  ...  in  (25).  By 
equating  powers  of  p  it  is  found  that 


2  2  1  (C1  '  ®b  )  {%  -  C2  )  4 

”l  '  P  -  8  A - 27  2  - 2\  ?  + 

®b  (C1  '  C2  ) 

V  V  Vl  c2 


(36) 


16 


Choosing  roots  to  correspond  to  outgoing  waves  or  bounded  functions  at 
infinity  gives 


Since  the  second  integrals  for  u  and  w  in  (23)  give  no  contribution 

until  t  >  —  c,  near  the  head  of  the  main  pulse  the  velocity  is  given 
2 

asymptotically  as  follows 


w 

T 


1 

2iri 


*^r 


,  P(t-c)  '  hd.p  c 
f  e  1  1 


(l  ♦  0(p2)) 


dp 


where  d1 


(38) 


By  the  convolution  theorem,  (38)  may  be  written 


where  Ai(-)  is  the  Airy  function  with  integral  representation  given  in 
Abromowitz  and  Stegun  ,  formula  10.4.32.  For  a  step  function  input, 
f(x)  =  h  (t) ,  w  is  simply  the  integral  of  the  Airy  function  and  is 
exponentially  small  for  5  >  t  and  oscillates  about  an  amplitude  of  1  for 
t  >5.  Fig.  2  shows  the  waveform  for  a  step  input.  Fig.  3  shows  the 
composite  waveform  schematically  for  both  the  leading  wave  and  the  main 
pulse.  Note  the  progressive  separation  and  decreasing  importance  of  the 
leading  fast  wave. 


18 

M.  Abromowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions, 
Dover,  New  York,  1965. 
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Figure  2:  Dispersion  in  Main  Pulse.  Linear  Case. 


Composite  Sketch  Showing  Separation  of  Leading  Wave  and  Main  Pulse.  Linear  Case 


The  initial  development  of  this  separation  was  shown  numerically  by 
Nunziato  and  Walsh  in  their  consideration  of  small  amplitude  waves  in 
porous  soils19  where  the  equations  are  formally  the  same  as  (20). 


V.  NONLINEAR  WAVE  EQUATIONS 


The  multiple  wave  structure,  which  the  linearized  equations  show  so 
clearly,  has  an  exact  counterpart  in  the  fully  nonlinear  theory.  The 
characteristic  wave  speeds  of  (5),  which  carry  discontinuities  in  w' ' 
and  w  or  u"  and  ii ,  are  given  by 


=  Sw,(w'  ,u,u‘) 

P2C22  =  Qu<(w'’u>u') 


(40) 


On  the  other  hand,  the  simplest  one-dimensional  theory  for  a  rod,  which 
is  given  by  (1) ,  and  which  may  be  obtained  from  (5)  by  setting  u'  =  0, 
F(w',u,o)  =  0,  has  characteristic  speeds  given  by 

■  *  <  V  • 

where  E  is  defined  in  (17) .  In  analogy  with  the  linear  theory  it  might 
be  expected  that  (41)  represents  a  low  order  approximation  for  wave 
speeds  in  the  main  pulse,  but  that  dispersive  effects  due  to  the  finite 
diameter  will  modify  the  pulse  shape.  Since  the  full  rod  theory,  (5), 
is  nonlinear,  it  is  to  be  expected  that  wave  speeds  will  be  modified  as 
well. 


Although  transform  methods  cannot  be  used  on  the  full  equations,  it 
is  possible,  by  suitable  scaling  and  choice  of  independent  variables,  to 
find  asymptotic  results  for  the  various  wave  orders  directly  from  the 
equations  themselves  rather  than  from  representations  for  exact  solutions. 
The  method  of  relatively  undistorted  waves,  introduced  by  Varley  and 
Cumberbatch^0  has  been  used  recently  by  Seymour  and  Parker10  to  examine 
the  leading  fast  wave  in  (5) .  The  same  approach,  also  called  the  method 
of  modulated  simple  waves,  may  be  used  to  investigate  the  structure  of 
the  main  pulse. 

To  simplify  the  analysis,  it  will  be  assumed  that  the  stored  energy 
in  (3)  has  the  decomposition 


E.  Varley  and  E.  Cumberbatch,  "Non-Linear  High  Frequency  Sound  Waves," 
J.  Inst,  Maths.  Applies.,  2,  1966,  pp.  133-143. 
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(42) 


W(w',u,u')  -  W^w’.u)  ♦  W2(u’) 


Let  the  dependent  variables  be  u,  v  «  w,  and  e  »  w* . 
are  rewritten  as 


Equations  (5) 


<5*  -  P  •  p , 


First  define  a  new  set  of  independent  variables 

C  «  6**  Z,  x  ■  6^  t,  (44) 

where  6  is  a  small  parameter.  This  scaling  implies  that  fixed  values  of 
C,  t  will  correspond  to  long  times  and  large  distances  as  6  tends  to  zero. 
Next  define  a  second  set  of  independent  variables  in  terms  of  c ,  x . 

a  -  a(C ,x) 

(45) 

z  =  6? 

In  terms  of  these  new  variables,  derivatives  with  respect  to  Z,  t  are 
replaced  as  follows. 

3  A  „  3  ,3/2  3 
3Z  6  C  3a  6  32 


3  A  3 
—  =  5  a  — — 

3 1  T  3a 

In  (45)  a  is  a  "fast"  and  z  a  "slow"  variable.  The  reason  for  the  odd 
combination  of  scaling  in  (44)  and  (45)  will  become  apparent  later.  It 
turns  out  that  a  has  properties  like  a  characteristic  variable.  Accord¬ 
ingly,  let  u>  and  k  be  defined  as 

w  •  aT  ,  <  ■  a^  .  (47) 

a  is  constant  along  curves  with  slope 


c  -  -  * 
dr  k 


It  is  convenient,  instead  of  u  and  k  to  use  the  slowness,  s  ■  —  ,  and 

1  c 

the  incremental  arrival  time  1  -  —  .  In  terms  of  these  variables,  the 
compatibility  condition  k t  becomes 
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V 


(49) 


s  -61 
a  z 


and,  taking  (42)  into  account,  equation  (43)  becomes 

sSuu,  *  clv»  *  *Sw'e*  •  St {Suuz  *  S«'ex} 


P  u  -  sP  ,v 
u  o  w'  a 


*  6&/-P  ,v  +  0  ,u  -  p0u  ) 
|  W*  Z  ^u'  2  TTTr 


(50) 


sv  +  e  *  6i,v 
a  a  z 


where  derivatives  in  t  and  c  are  retained  on  the  right  hand  side  of  (50)2 
for  compactness. 


To  lowest  order  the  right  hand  sides  of  (50)  may  be  set  equal  to 
zero.  This  requires  that  the  determinant  of  coefficients  of  ua,  vq,  e 
must  vanish  giving  ‘ 


1 

7 


(SI) 


and  u  ,  v  ,  e  must  occur  in  the  ratio 
a  a  a 

V  Va:eo  =  -  v:  -  V  1  •  (52) 

With  the  energy  partition  of  (42)  neither  S  nor  P  depends  on  u'.  Thus 
(52)  may  be  integrated  to  find  two  of  u,v,e  in  terms  of  the  third.  For 
example 


L 


V  ■  v(u) 
c  ■  e(u) 


(53) 


The  right  hand  side  of  (50)  must  satisfy  a  compatibility  condition.  Thus 
multiplication  of  (50)  by  the  left  proper  vector 


(-1,  sv,  sSw,) 


(54) 


gives  the  transport  equation 


_  „  \  21  9  [  s  3/s 

Sz  -  sEvz  "  ^  I c2  liJ-  Lila  (l 

-  r  if  tr  dKi  Mil' 

2  1 

where  P2c2  •  *»  0  , .  Since  S,s,E  and  v  all  depend  on  u  through  (53),  the 

left  hand  side  of  (55)  may  be  written  as  #(u)uz  so  that  the  dependent 

variable  is  u  alone.  However,  the  mapping  a(c«t)  must  be  developed 
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“*)] 


(55) 


simultaneously  with  the  solution  u(a,z)  since  the  incremental  arrival 
time  i  is  also  unknown. 

The  properties  of  (55)  will  not  be  studied  further  here  except  for 
the  following  two  remarks.  First,  the  linearized  version  of  (55)  can  be 
shown  to  have  solutions  like  (39) ,  and  second,  for  small  (but  finite) 
amplitude  waves,  the  transport  equation  reduces  to  the  Korteweg-de  Vries 
equation.  This  latter  fact  may  also  be  deduced  directly  from  (5)  follow¬ 
ing  the  prescription  of  Leibovich  and  Seebass  .  Many  other  authors  have 
also  noted  that  the  Korteweg-de  Vries  equation  describes  longitudinal 
waves  in  thin  rods,  eg.,  Nariboli  and  Sedov.22 


VI.  STEADY  WAVES 


Equation  (49)  indicates  that  if  £  =  0(1),  then  s,  and  hence  c^,  is 

almost  constant.  Therefore,  some  solutions  of  (55)  should  be  nearly 
steady  waves,  perhaps  approaching  them  asymptotically  in  time.  Accord¬ 
ingly,  the  properties  of  steady  solutions  to  (5)  will  be  examined. 

A  steady  wave  is  one  in  which  the  solution  to  (5)  may  be  written 


w  =  w(Z  -  ct) 
u  »  u(Z  -  ct) 


(56) 


where  c  is  an  unspecified,  but  constant,  wave  speed, 
become 


S 


w" 


» 


Q'  -  P 


P2c 


u’ 


The  equations 


(57) 


where  the  prime  now  denotes  differentiation  with  respect  to  Z  -  ct. 


2 1 

S.  Leibovich  and  A.  R.  Seebass,  "Examples  of  Dissipative  and  Dispersive 
Systems  Leading  to  the  Burgers  and  the  Korteweg-de  Vries  Equations," 
in  Nonlinear  Waves,  S.  Leibovich  and  A.  R.  Seebass,  eds.,  Cornell 
University  Press,  Ithaca,  N.Y. ,  1974. 

22 

G.  A.  Nariboli  and  A.  Sedov,  "Burgers' s-Korteweg-de  Vries  Equation  for 
Visscoelastic  Rods  and  Plates,"  J.  Math.  Anal,  and  Appl.,  32,  1970, 
pp.  661-667. 
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There  are  two  integrals  of  the  motion. 


S  ■  PjC  w'  ♦  A 


Q  u'  ♦  A  w*  -  W 

-I'2(»2“'2  -  »l“'2)  * 


(58) 


where  A  and  B  are  constants.  Thd  second  integral  is  found  by  multiplying 

(57).  by  w'  and  (52),  by  u' ,  adding  the  two  equations,  and  noting  that 

^  /  3W  aw  \  i  I  2/  2  2) 1 

the  result  may  be  written  as  (  g^jr  u'  +  g^r  w'  -  W)  *  c^jW’  -^u'  /  . 

After  integration  (58) ^  is  substituted  into  this  expression  to  give  (58)  2 


If  (58)  can  be  solved  for  w*  and  u* ,  they  may  be  expressed  as 
follows. 


w'  =  F(u;  A,B,c) 
u'2  =  G(u;  A,B,c) 


(59) 


2 

Equation  (59),  is  written  in  terms  of  u'  rather  than  u'  since  both  of 

(58)  are  even  in  u' .  Now  (59) 2  may  be  solved  by  quadrature  and  then 

(59) .  gives  the  distribution  in  longitudinal  strain  parametrically 
through  u. 


To  examine  the  qualitative  nature  of  solutions  and  the  role  of  the 
constants  A,  B,  and  c,  first  consider  the  case  A  =  B  =  0.  For  ease  of 
analysis  the  strain  energy  will  be  assumed  to  be  partitioned  as  in  (42) . 
For  small  but  finite  amplitude  waves  let  W  be  expressed  as  a  power  series 
First  note  that 


where  W  is  the  three-dimensional  strain  energy  per  unit  volume.  With 
(12)  used  in  (60) ,  W  may  be  expressed  as 

W  =  j  (A+2p)  w'2  +  2Awuw'  +  2(X+y)u2 

+  j  a2  yu'2  +  CjU3  ♦  C2u2w'  (61) 

+  C3uw»2  +  C4w'3  +  _  , 
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where  the  C's  are  combinations  of  higher  order  elastic  moduli. 
Equation  (58)  ^  with  A  =  0  gives 


w'  = - Q - -  u  +  0(u2)  , 

>.  +  2 y  -  pc 

and  (58) 2  with  A  =  B  =  0  and  the  use  of  (62)  gives 

2  2  3 

u'  =  Cu  +  Du  ♦  ...  , 

where 

c  =  8(A+u)  (E-pc2) 

a2(y-pc2)(X+2u-pc2) 


(62) 


(63) 


(64) 


and  D  may  be  either  positive  or  negative,  depending  on  the  exact  combina¬ 
tion  of  higher  order  elastic  moduli.  C  may  be  either  positive  or  negative 
as  well. 


C  >  0  if  cfc  <  c  <  cx 
C  <  0  if  c2  <  c  <  c^ 


(65) 


wheTe  Cj,  c2,  and  are  the  linear  elastic  wave  speeds  in  this  context. 

2 

Figure  4  shows  sketches  of  u'  vs.  u  and  u'  vs.  u  for  the  various 
combinations  of  signs  of  C  and  D.  The  closed  loops  in  (4e)  and  (4f) 
correspond  to  solitary  waves  that  are  symmetrical  in  Z  about  the  point  of 
maximum  amplitude  and  extend  from  Z  *  -®  to  Z  *  +®  .  The  open  curves  in 
(4e)  -  (4h)  correspond  to  waves  in  which  the  amplitude  grows  without 
limit.  Clearly,  to  be  useful,  solutions  of  this  kind  must  join  onto 
unsteady  solutions  before  the  amplitudes  become  too  extreme. 


Now  consider  the  case  where  A  =  0  but  B  +  0.  Instead  of  (63)  the 
following  holds. 

u'2  =  B  +  Cu2  +  Du^  (66) 

The  constant  B  is  proportional  to  B,  but  has  opposite  sign.  B  merely 
shifts  the  curves  in  (4a) -(4d)  vertically.  Figure  5  shows  two  of  the 
possibilities.  The  closed  loops  in  (5c)  and  (5d)  correspond  to  periodic 
solutions. 


Next  consider  the  case  when  both  A  i  0  and  B  /  0.  The  integrals  of 
the  motion  may  be  rewritten  as  follows 
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ase  Plane  Sketches  for  Steady  Waves,  A=B=0  in  every  case, 
bl  and  F>  C>0.  D<0.  In  c)  and  g)  C<0,  D>0.  In  d)  and  h) 


B>0) ,  otherwise  like  4c). 


(67) 


S  -  SQ  -  Plc2(w'.w^) 

Qu'  ♦  S0(W-wp  -  (W  -  WQ) 

-  \  c2{p2u*2  -  Pl(w'-w;)2^  B  , 

where  Wq  =  W(w',u  ,0),  etc.  These  integrals  are  the  same  as  before,  but 
now  they  are  centered  on  the  point  (w^,uq,0)  and  SQ  plays  the  role  of  A. 
The  character  of  the  solutions  is  exactly  the  same  as  in  Figures  4  and  5, 
but  now  the  curves  are  shifted  to  the  right  (SQ<0)  or  left  (SQ>0) . 

Finally,  suppose  that  the  elastic  constants  are  such  that  D  «  0. 
Instead  of  (66) ,  equation  (68)  must  be  used. 

u'2  «  B  +  Cu2  +  Hu4  (68) 


This  should  occur  when  the  expansion  is  centered  at  an  inflexion  point  of 
the  static  stress-strain  curve.  C  is  the  same  as  before  and  H  is  deter¬ 
mined  by  higher  order  elastic  moduli.  Some  of  the  possibilities  are 
shown  in  Figure  6.  Either  positive  or  negative  solitary  waves  are 
possible  according  to  (6c) ,  and  (6d)  is  a  periodic  solution  of  unusual 
shape . 

The  dynamic  trajectories,  D,  of  the  nonlinear  steady  waves  may 
deviate  considerably  from  the  static  curve  P(w',u,0)  =  0  in  the  w'-u 
plane.  To  see  this,  again  suppose  that  the  strain  energy  is  partitioned 
as  in  (42)  and  consider  the  case  A  *  0.  Figure  (7a)  shows  a  contour 
plot  of  on  the  w'-u  plane  with  the  curve  P  =  =0  and  several  lines 

of  u  =  const  drawn  in.  In  Figure  (7b)  several  curves  of  S  *  W,  vs.  w' 

Aw' 

with  u  =  constant  are  shown.  Shown  as  a  dashed  curve  on  (7b)  is  the 
static  stress-strain  curve,  which  has  somewhat  smaller  slope  than  any  of 
the  curves  for  u  =  const,  according  to  (17).  From  reference  to  (58) ^  it 
can  be  observed  that  points  of  the  nonlinear  trajectories  lie  on  the1 
intersection  of  the  surfaces  =  S(w',u)  and  S2  =  pjC2w'.  In  the 

projection  of  Figure  (7b)  S2  is  a  straight  line^through  the  origin  with 
slope  P^c  ,  for  example  the^line  p^jp2  (for  c  less  than  Cv,  the  bar  speed 
at  w*  =  0) .  By  comparing  the  line  pjP2  curves  For  S  at  u  *  const 

and  the  static  stress/strain  curve,  it  can  be  seen  that  between  points 
Pj  and  p2  the  dynamic  trajectory  deviates  to  values  of  u  less  than  those 
on  P  =  o,  and  at  points  outside  the  segment  p.p2  the  trajectory  deviates 
to  greater  values  of  u  than  those  on  P  =  0.  thus,  the  trajectory,  D,  is 
qualitatively  like  the  dashed  curve  in  Figure  (7c) . 

Maximum  or  minimum  amplitudes  of  u  will  occur  where  u'  =  0.  With 
u'  =  0,  equation  (58) 2  generates  curves,  C,  in  the  w'-u  plane  that 
correspond  to  the  intersection  of  the  surfaces  =  W^(w',u)  and 
W2  =  j  P^c2w'  -B.  The  intersection  of  curves  C  and  D  mark  the  ends  of 
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ane 


Figure  7:  Comparison  of  Steady  Dynamic  Trajectory  with  Static  Trajectory 
a)  Contours  of  *  const,  and  Static  Curve  P=0.  b)  Intersec¬ 
tion  of  Surfaces  S(w',u)  and  p.c^w' .  c)  Dynamic  Trajectory 
Deduced  from  4b).  d)  Intersection  of  Cylinder  l/2p,c^w'^-B 
with  Widest  Projection  of  W7(w',u). 


the  dynamic  trajectory.  By  differentiating  (58)2  with  respect  to  w',  it 
may  be  seen  that  the  points  on  C  satisfy  ^ 


P 


S 


2  . 
PjC  w' 


(69) 


At  the  intersection  of  C  and  D,  (58) ^  holds  as  well  so  that  at  those 
points  1 


0 


(70) 


That  is,  either  P  =  0  or  the  slope  of  C  is  horizontal. 

Now  consider  further  the  intersection  of  surfaces  and  W2  with 
B  <  o  (B  >  o) .  Fig  (7d)  shows  the  projection  parallel  to  the  u  axis  of 
the  widest  points  of  =  W^(w',u).  These  points  correspond  to  the 
curve  P  =  0  in  (7a),  and,  in  fact,  the  projected  potential  well  in  (7d) 

-  1  2  2 

generates  the  static  stress/strain  curve.  The  surface  J  plc  w'  “  ® 

is  a  parabolic  cylinder  with  generators  parallel  to  the  u  axis.  An  end 
view  of  the  cylinder  is  shown  as  the  dashed  parabola  in  (7d) .  Since  the 
level  curves  of  are  closed,  as  shown  in  (7a),  and  since  the  level 
curves  of  W2  are  straight  lines  parallel  to  the  u-axis,  the  curves  C  are 
qualitatively  similar  to  those  sketched  in  (7c).  Thus,  for  the  case 
considered,  which  corresponds  to  Figure  (Sd) ,  the  motion  lies  along  D 
and  either  oscillates  between  points  p3  and  p4  or  begins  at  p,  and  moves 
to  the  right  until  interrupted  by  an  unsteady  motion.  In  either  case, 
note  the  divergence  of  D  from  the  curve  P  *=  0. 

Other  cases  for  different  shapes  of  the  potential  well  W1(w’,u)  may 
be  considered  in  a  similar  manner.  However,  from  teh  preceding  qualita¬ 
tive  discussion  it  is  clear  that  te  cubic  expansion  in  equation  (63)  leads 
to  the  qualitatively  correct  curves  in  the  phase  plane  as  shown  in 
Figures  (4)  and  (5),  provided  that  the  function 

K(u'2)  =  Q(u')u'  -  W2(u')  -  \  p2c2  u'2 

2 

is  monotonic  and  hence  single  valued  in  u'  .  If  this  is  not  the  case, 
then  the  curves  of  u'^  vs.  u  in  Figures  (4),  (5),  and  (6)  will  have 
bifurcation  points.  The  analysis  would  then  be  somewhat  more  complicated 
and  will  not  be  attempted  here. 


VII.  REMARKS 

This  report  has  been  motivated  by  the  need  to  interpret  wave  propa¬ 
gation  experiments  in  rods  and  bars.  These  experiments  are  of  two  types. 
In  one,  a  cylindrical  bar  is  fired  from  an  air  gun  so  as  to  impact  an 
identical  bar,  which  is  at  rest  in  an  evacuated  chamber.  The  bar  at  rest 
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is  instrumented  at  several  stations  along  its  length  with  foil  resistance 
strain  gages.  In  the  other  experiment,  a  long  rod  penetrator,  which  is 
also  instrumented  with  strain  gages,  is  impacted  by  a  target  plate.  3 
In  either  case  a  transient  plastic  wave  propagates  into  the  bar,  strain¬ 
time  data  are  recorded  for  the  instrumented  stations  of  the  bar,  and 
then  the  data  are  interpreted  with  the  aid  of  equations  (2) .  By  compar¬ 
ing  with  streak  camera  records  Hauver  has  shown  that  (2^  gives  very 
accurate  results  for  particle  velocity  ,  but  (2)2  always  gives  higher 
stresses  than  the  quasi-static  stress-strain  curve  even  though  the  high 
strength  steels  and  aluminums  used  are  not  rate  sensitive  in  Hopkinson 
bar  tests.  4  The  only  other  source  of  excess  stress  that  has  been 
suggested  is  the  effect  of  finite  transverse  dimensions,  hence  the 
present  interest  in  a  nonlinear  theory  that  accounts  for  radial  motion. 

This  report  has  dealt  only  with  nonlinear  elastic  materials,  not 
elastic-plastic  materials.  Nevertheless,  there  are  several  results  that 
are  at  least  suggestive  of  general  features  that  should  apply  to  the 
elastic-plastic  case.  These  are  as  follows. 

a.  The  averaging  in  equation  (10)  and  the  identifications 
made  in  equations  (11)  and  (12)  are  independent  of  constitutive  rela¬ 
tions.  Therefore,  the  equations  of  motion  (5)  hold  for  the  elastic- 
plastic  case  as  well  as  the  elastic  case,  but  of  course  equations  (6) 
will  not  be  appropriate  then. 

b.  With  v  ■  w  and  e  ■  w*  equations  (1)  are  entirely  consistent 
with  (5) .  As  a  consequence,  if  each  level  of  strain  propagates  at  a 
constant  velocity,  experimental  data  may  still  be  interpreted  through 
the  use  of  (2)  for  both  particle  velocity  and  stress.  In  this  integra¬ 
tion  no  use  is  made  of  (5)2  so  that  radial  strain  u  remains  undetermined. 
In  particular,  there  is  no  requirement  for  u  to  have  equilibrium  values 
corresponding  to  the  strain  e.  Thus,  the  computed  stress  S  should  be  as 
accurate  as  the  computed  particle  velocity  w  even  though  S  cannot  be 
measured  independently  as  can  w.  In  a  compressive  wave  the  excess 
stress  may  be  accounted  for  by  underexpansion  in  the  radial  direction 
caused  by  inertial  lag. 

c.  In  impact  experiments  of  the  type  described  here  it  is 
commonly  observed  that  the  plastic  wave  speed  becomes  nearly  constant 
at  the  higher  values  of  strain.  Thus,  steady  waves  are  of  interest  at 
least  as  asymptotic  solutions.  Because  of  dissipation  it  is  probably 
unreasonable  to  expect  that  solitary  waves  or  periodic  waves  will  occur 
in  elastic-plastic  materials,  but  pure  loading  waves  should  be  similar 
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in  either  elastic  or  elastic-plastic  materials.  Thus,  unbounded  waves 
of  the  type  shown  in  Figures  (4g),  (4h) ,  (5c),  or  (5d)  could  occur.  In 
an  actual  physical  problem  such  a  wave  must  be  interrupted  by  an 
unsteady  wave  or  perhaps  a  free  boundary  (as  in  the  penetration  experi¬ 
ment)  . 

Future  work  will  attempt  to  elaborate  on  some  of  these  issues. 
Initial-boundary  value  problems  will  be  posed  with  an  emphasis  on 
asymptotic  solutions.  The  elastic-plastic  case  will  also  be  formulated 
and  the  characteristic  features  of  those  solutions  examined  and  compared 
with  elastic  solutions. 
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